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We have investigated low-temperature transport properties of two-dimensional arrays of 
superconductor-normal-metal-superconductor (SNS) junctions. It has been found that in two- 
dimensional arrays of SNS junctions (i) a change in the energy spectrum within an interval of the 
order of the Thouless energy is observed even when the thermal broadening far exceeds the Thouless 
energy for a single SNS junction; (ii) the manifestation of the subharmonic energy gap structure 
(SGS) with high harmonic numbers is possible even if the energy relaxation length is smaller than 
that required for the realization of a multiple Andreev reflection in a single SNS junction. These 
results point to the synchronization of a great number of SNS junctions. A mechanism of the 
SGS origin in two-dimensional arrays of SNS junctions, involving the processes of conventional and 
crossed Andreev reflection, is proposed. 
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Andreev reflection is a microscopic mechanism respon- 
sible for the charge transport through a normal-metal- 
superconductor interface Qj. An electron- like quasipar- 
ticle in the normal metal with an energy lower than the 
superconducting gap (A) is reflected from the bound- 
ary as a hole-like quasiparticle, while a Cooper pair is 
transferred to the superconductor. If the normal metal 
is sandwiched between two superconductors, there is an 
additional charge transfer mechanism, namely, multiple 
Andreev reflection (MAR). This mechanism was pro- 
posed in [2j to explain the subharmonic energy gap struc- 
ture (SGS) observed as dips in the differential resistance 
curves at voltages eV n = 2A/n (n is an integer). The 
MAR concept is as follows: owing to sequential Andreev 
reflections from normal-metal-superconductor (NS) in- 
terfaces, a quasiparticle passing through the normal re- 
gion can accumulate an energy of 2 A, which is suffi- 
cient for the transition to single-particle states of the 
superconductor. Although, by now, ample experimental 
data on the transport properties of NS and SNS junc- 
tions are available [3j and the existing theoretical models 
fairly well describe the phenomena observed in the ex- 
periments 0, 0, 01 1 the properties of multiply connected 
SNS systems are poorly investigated. It should be noted 
that systems consisting of superconducting islands incor- 
porated into a normal metal are spontaneously formed 
in disordered superconducting films • In view of these 
circumstances, it is of interest to study the properties of 
model multiply connected SNS systems with supercon- 
ducting and normal regions formed in a controlled way. 

This paper is devoted to the properties of two- 
dimensional arrays of SNS junctions fabricated on the 
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basis of a 20-nm-thick PtSi superconducting film (with a 
critical temperature of T c = 0.64 K) The transport 
parameters of the initial PtSi film were as follows: the re- 
sistance per square at T = 4.2 K was R sq = 22.8 fi, the 
mean free path was I = 1.35 nm, and the diffusion coeffi- 
cient was D = 7.2 cm 2 /s. The initial samples were fabri- 
cated by photolithography in the form of Hall bars 50 /im 
wide and 100 /mi long. Then, by electron beam lithogra- 
phy with subsequent plasma chemical etching, the initial 
film was thinned in preset regions. Figure schemat- 
ically represents the structure under study. It consists 
of periodically arranged film areas with a thickness of 
20 nm (islands), between which the film is thinned by 
plasma chemical etching. The period of the structure 
is 1 /im, and the island dimensions are 0.8 x 0.8 /tm. 
This structure completely covers the whole area of the 
Hall bar (Fig. [TJ>). Thus, the number of islands be- 
tween the potential terminals is 50 x 100. To control 
the depth of etching, a reference film was etched uni- 
formly over its surface simultaneously with the structure 
etching. The low-temperature transport measurements 
were performed by a standard four-probe low-frequency 
(~10 Hz) technique. An ac current was within the range 
1 - 10 nA. 

As the experiment shows, the decrease of film thick- 
ness by plasma chemical etching leads to the suppression 
of the critical temperature. At the same time, in the is- 
lands, whose thickness remains equal to that of the initial 
PtSi film, the suppression of T c is minor. Hence, there 
is a certain temperature region within which the thinned 
film is in the normal (N) state and the islands are in the 
superconducting (S) state. In this temperature region, 
the structures under study represent a two-dimensional 
array of SNS junctions. Moreover, since the supercon- 
ducting and normal regions of the SNS junctions, fab- 
ricated in the aforementioned way, consist of the same 
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FIG. 1: Schematic representation of the structure: (a) profile 
of the film modulated in thickness (the geometric dimensions 
are indicated in the figure); (b) position of the structured film 
on the sample used for the measurements. 

material, the formation of tunnel barriers at the NS in- 
terfaces is excluded and a high transparency of the NS 
interfaces can be a priori expected. 

Figure [21 shows the temperature dependence of resis- 
tance for two samples that only differ in the film thick- 
ness between the islands. The resistance per square of 
the reference film at T = 4.2 K is 995 il for sample S#l 
and 1483 f2 for sample S#2. For all structures studied, 
a noticeable decrease in resistance with decreasing tem- 
perature is observed slightly below the temperature T c of 
the initial 20-nm-thick PtSi film. As the temperature de- 
creases (Fig. [3J|, a small decrease in the resistance is first 
observed at T ~ 0.61 K; then, the resistance increases 
reaching a maximum and then decreases again. Such an 
anomalous behavior of the temperature dependence of 
resistance in the narrow temperature range near T c was 
never observed before, and the origin of this behavior is 
unknown. 

To explore the nature of this anomaly, the current- 
voltage characteristics of the samples have been studied. 
Figure shows the dependences of the differential resis- 
tance (dV/dl) on the bias voltage for the sample S#l. 
The measured values of dV/dl are given per square. The 
total voltage is divided by the number of junctions (100) 
in the rows between the potential probes. For a square 
array, this procedure yields the average voltage across 
one SNS junction, and this quantity is represented by the 
abscissa axis in Fig. (below, all dependences of the dif- 
ferential resistance on the bias voltage are plotted in the 
same coordinate systems) . The dependences of dV/ dl-V 
are symmetric with respect to the direction of the cur- 
rent. As one can see from the dependences shown in 
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FIG. 2: Temperature dependence of the resistance per square 
for samples S#l and S#2. The dash in the plot indicates the 
superconducting transition temperature for the initial 20-nm- 
thick PtSi film. The arrows indicate the temperature range 
from 0.545 to 0.599 K. 

Fig. [3 an increase in temperature leads to a suppression 
of the excess conductivity and, at T ~ 570 mK, the min- 
imum observed at V = is replaced by a maximum. It 
should be emphasized that the change of a minimum to 
a maximum in the differential resistance dependences at 
zero bias voltage occurs in the same temperature region 
within which the anomaly is observed in the temperature 
dependence of resistance (in Fig. [2 for sample S#l, this 
interval is indicated by arrows). 

Let us estimate the characteristic energy scales for the 
two-dimensional array of SNS junctions under study. The 
energy corresponding to the voltage, at which the min- 
imum is suppressed at zero bias and temperatures be- 
low 570 mK, is estimated as eV ex ~ 3 ^eV. The pres- 
ence of the maximum in the differential resistance at 
T > 570 mK may be caused by the minimum that occurs 
in the density of states in the normal region. The theo- 
retical study of the properties of diffusive NS junctions 
shows that the decrease in the density of states is caused 
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FIG. 3: Differential resistance per square versus the bias volt- 
age per SNS junction for sample S#l. All dependences, ex- 
cept for the lower one, are sequentially shifted upwards by 
0.5fi. 



by the presence of coherent Cooper pairs in the normal 
metal. According to the theory, at a distance x from 
the NS interface, only the pairs from the energy window 
Ef ± £rh (the coherence window), where £rh = hD/x 2 
is the Thouless energy, will remain coherent. Therefore, 
the density of states has a minimum at energies E < E^ 
and a maximum at E ~ i/rh, and value of the density of 
states in the maximum is higher than that in the normal 
metal. The experimental study of the density of states as 
a function of the distance to the NS interface ^} showed 
a good agreement with the theory. Note that this ef- 
fect manifests itself in two ways. On the one hand, the 
presence of coherent Cooper pairs leads to a decrease 
in the resistance of the normal region (an analog of the 
Maki-Thompson correction), and, on the other hand, the 
pairing leads to a decrease in the density of single-particle 
states (an analog of the correction to the density of states 
in the Cooper channel). These competing contributions 
may lead to a nonmonotonic temperature dependence of 
the resistance of SNS junctions Let us estimate the 
Thouless energy Et\i — HD/L 2 for the structure under 
study, where L = 0.2 /jm is the length of the normal re- 
gion and D = 2 cm 2 /s is the diffusion coefficient for the 
reference film. Then, the Thouless energy is i?Th — 3 
fJ.eV, which nearly coincides with the characteristic volt- 



age value (in Fig. this value is indicated by the ver- 
tical straight line). However, for the observation of the 
effects associated with changes in the quasiparticle spec- 
trum at the Thouless energy, the necessary condition is 
kT < i?Th- In the temperature range of interest, e.g., at 
T = 580 mK, the thermal broadening is kT ~ 50 /ieV, 
which noticeably exceeds the correlation energy Erh for 
the SNS junctions forming the array. Thus, one should 
expect that the processes determined by E^ would not 
manifest themselves in the differential resistance curves. 
However, the experimental dependences clearly display 
the suppression of the maximum at a voltage of about 
£>Th/e. This fact suggests that the effect observed in the 
experiment is collective; i.e., the whole array is responsi- 
ble for its manifestation. Note that the total bias voltage 
is eV ex ■ 100 = eVs = 300 fieV and it is much higher 
than kT. We will return to discussing this issue after 
presenting results that, in our opinion, also testify to a 
correlation in the behavior of the array of SNS junctions. 

Figure EJi shows the dependences of the differential re- 
sistance on the bias voltage per one SNS junction. These 
dependences have pronounced minima at some voltage 
values (the minima are marked by a, /3, and 7 in the 
plot). When temperature increases, the minima are 
shifted to lower voltages and, in the interval from 550 
to 600 mK, one more minimum (5) appears; at lower 
temperatures, this minimum is absent. The tempera- 
ture dependences of the positions of these features are 
shown in Fig. 0Jd by different symbols. The solid lines in 
Fig. 2t> represent the temperature dependences of the su- 
perconducting energy gap 2A(T)/en, where A(T) is the 
temperature dependence of the superconducting energy 
gap predicted by the BCS microscopic superconductiv- 
ity theory and n is an integer (the corresponding values 
are indicated in Fig. 0Jd). The exact value of A(0) for 
platinum monosilicide is unknown. The estimate by the 
formula A(0) = 1.76fcT c (T c = 0.88 K for bulk PtSi [H) 
yields the value A/e ~ 133 fxV. If we set A/e = 126 /xV, 
the positions of the features in the voltage dependences of 
the differential resistance will be determined by the con- 
dition eV = 2A(T)/n. From Fig.QJs, one can see that the 
temperature dependences of the positions of minima ob- 
served in the experimental differential resistance curves 
functionally coincide with the dependences 2A(T)/en. 
From this fact, we can conclude that these features rep- 
resent the subharmonic energy gap structure (SGS). 

However, there are two circumstances clearly indicat- 
ing that it is difficult to treat the features observed at 
bias voltages equal to 2A(T) as the manifestation of the 
SGS caused by multiple Andreev reflections that occur 
independently in each SNS junction and that this inter- 
pretation is impossible for large values of n. First, the 
conventional mechanism of the SGS formation in a sin- 
gle SNS junction implies that, for the appearance of a 
subharmonic of number n, an n-fold passage of quasipar- 
ticles through the normal region is necessary without any 
energy relaxation. In other words, the energy relaxation 
length (l e ) must be no smaller than n ■ L 6j. The en- 
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FIG. 4: Sample S#l: (a) differential resistance per square 
versus the bias voltage per SNS junction (i.e., the voltage 
across the potential terminals is divided by 100) ; (b) positions 
of the features marked in Fig.^K as functions of temperature. 
The solid lines represent 2A(T)/en. The respective values of 
n are indicated in the plot. The dependences are calculated 
for A(0) = 126 fieV and T c = 0.625 K. 



ergy relaxation length can be estimated from the known 
phase breaking length l v . In the temperature range of 
interest, l e can exceed l v by no more than an order of 
magnitude ^3|. When a carrier acquires an energy of 
2A = 252 fieV, it is heated to a temperature of ~3 K. Ac- 
cording to our magnetotransport measurements, at this 
temperature, we have l v — 40 nm. Then, the most op- 
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FIG. 5: Andreev reflection process in a semiconductor repre- 
sentation for the structure consisting of alternating supercon- 
ducting and normal regions. 



timistic estimate yields l t ~ 0.4 nm. This value is much 
smaller than 14 x 0.2 //m = 2.8 /im, which is neces- 
sary for the realization of the subharmonic with n = 14. 
Second, it is difficult to expect that normal regions are 
fully identical, i.e., characterized by exactly the same re- 
sistance. Evidently, we deal with a network of different 
resistances. Therefore, when a current flows through the 
structure, the voltages across the normal regions are dif- 
ferent. This means that, if, for some of the normal re- 
gions, the condition eV — 2A/n is satisfied at a given 
value of the current, for other regions it will be satisfied 
at some other values of the current. Only when the resis- 
tances are close to each other, it is possible to observe the 
subharmonic energy gap structure as a result of the sta- 
tistical averaging; however, this structure will be smeared 
in proportion to the deviations of actual resistances from 
the average value. These speculations lead to a natural 
conclusion that, for a network of random resistances, the 
observation of the subharmonic energy gap structure is 
impossible. Contrary to this conclusion, the experiment 
shows that the SGS still manifests itself, although it is 
somewhat irregular; i.e., it lacks some of the harmon- 
ics. Similar results for an array of SNS junctions of other 
configuration were obtained in |14| . By now, no theory 
has been developed to describe multiply connected SNS 
systems, and, on the basis of the theoretical results ob- 
tained for single SNS junctions, it is difficult to explain 
the observation of subharmonics with numbers reaching 
n = 14. 

In considering the properties of single SNS junctions, 
the superconducting regions are usually assumed to be 
large and serving as reservoirs for electrons. The funda- 
mental distinction of the array of SNS junctions is the 
finite size of the superconducting regions (is)- There- 
fore, it is necessary to take into account the possibility 
of the mutual influence of the Andreev reflection pro- 
cesses that occur at the interfaces of one superconduct- 
ing island. In NSN junctions, in addition to the "com- 
mon" Andreev reflection, when an electron-like excita- 
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tion is reflected at the interface of the superconductor 
with transformation into a hole-like excitation (in Fig. [5] 
process 1: ej hi), an additional process of crossed An- 
dreev reflection a <-> ft^ (process 2) is possible, when 
the electron-like and hole-like excitations are on different 
sides of the superconducting region 0, 0, 0. As a 
consequence, such a reflection in combination with the 
common Andreev reflection, which occurs on one side of 
the superconductor (process 3), leads to the passage of 
quasiparticles without any energy loss. Note that, for the 
observation of the nth subharmonic of Andreev reflection, 
the realization of this mechanism of charge transport does 
not require that the voltage across each of the normal re- 
gions be exactly equal to V n = 2 A/en. It is sufficient 
that the total voltage across the chain of normal regions 
connected in series through superconducting regions be 
equal to 2A/e: 

i+n—l 

J2 V s = 2A/e. 

s—i 

Moreover, for the observation of the nth subharmonic 
of Andreev reflection, the condition l e > nL needs not 
necessarily be satisfied. It is sufficient that the energy 
relaxation length be greater than the distance between 



the superconducting regions (l e > L). We also note that 
the crossed Andreev reflection process makes a consid- 
erable contribution to the charge transport through an 
NSN junction even if the size of the superconducting re- 
gion is several tim es g reater than the superconducting 
coherence length £ (TEl Il8j. In our case, £(0) ~ 70 nm, 
and at T = 100 mK, we have L s /£ ~ 10. 

We believe that this charge transport mechanism is 
responsible for the features observed in the experiment. 
Returning to the behavior of the two-dimensional array 
of SNS junctions near zero bias, we note that the syn- 
chronous Andreev reflection on both sides of the super- 
conducting regions should lead to correlated changes of 
the energy spectra in normal regions by analogy with the 
formation of minibands and minigaps in semiconductor 
superlattices. 

We are grateful to M.V. Feigel'man for useful discus- 
sions and for the interest taken in our work. This work 
was supported by the Russian Academy of Sciences (the 
program "Quantum Macrophysics" ) , the Ministry of Sci- 
ence of the Russian Federation (the program "Super- 
conductivity of Mesoscopic and Strongly Correlated Sys- 
tems"), and the Russian Foundation for Basic Research 
(project no. 03-02-16368). 



[1] A.F. Andreev, Zh. Eksp. Teor. Fiz. 46, 1823 (1964) [Sov. 
Phys. JETP 19, 1228 (1964)]. 

[2] T.M. Klapwijk, G.E. Blonder, and M. Tinknam, Physica 
109-110B+C, 1657 (1982). 

[3] B. Pannetier and H. Courtois, J. of Low Temp. Phys. 
118, 599 (2000). 

[4] G.E. Blonder, M. Tinkham, and T.M. Klapwijk, Phys. 
Rev. B25, 4515 (1983). 

[5] M. Octavio, M. Tinkham, G.E. Blonder, and T.M. Klap- 
wijk, Phys. Rev. B27, 6739 (1983); K. Flensberg, 
J. Bindslev Hansen, M. Octavio, ibid. 38, 8707 (1988). 

[6] E.V. Bezuglyi, E.N. Bratus', V.S. Shumeiko, 
G. Wendin, and H. Takayanagi, Phys. Rev. B62, 
14439 (2000). 

[7] M. V. Feigel'man, A. I. Larkin, and M. A. Skvortsov, 

Phys. Rev. Lett. 86, 1869 (2001). 
[8] R. A. Donaton, S. Jin, H. Bender, et al, Microelectronic 

Engineering 37/38, 504 (1997). 
[9] B.J. van Wees, P. de Vries, P. H. C. Magnee, and 

T. M. Klapwijk, Phys. Rev. Lett. 69, 510 (1992). 



[10] S. Gueron, Norman O. Birge, D. Esteve, and M. H. De- 

voret, Phys. Rev. Lett. 77, 3025 (1996). 
[11] Anatoly F. Volkov and Hideaki Takayanagi. Phys. Rev. 

B 56, 11184 (1997). 
[12] B.T.Matthias, T. H. Geballe, and V.B.Compton. Rev. 

Mod. Phys. 35, 1 (1963). 
[13] S.I. Dorozhkin, F. Lell and W. Schoepe, Solid State 

Comm. 60, 245 (1986). 
[14] T.I. Baturina, Z.D. Kvon, and A. E. Plotnikov, Phys. 

Rev. B 63, 180503(R) (2001). 
[15] J.M. Byers and M.E. Flatte, Phys. Rev. Lett. 74, 306 

(1995). 

[16] G. Deutscher and D.Feinberg, Appl. Phys. Lett. 76, 487 
(2000). 

[17] Arne Jacobs and Reiner Kiimmel, Phys. Rev. B 64, 
104515 (2001). 

[18] D. Beckman, H. B. Weber, H. v. Lohneysen, Phys. Rev. 
Lett. 93, 197003 (2004). 



